!===============================================================================
! OLAM version 3.3

! Copyright (C) 2002-2008; All Rights Reserved; 
! Duke University, Durham, North Carolina, USA 

! Portions of this software are copied or derived from the RAMS software
! package.  The following copyright notice pertains to RAMS and its derivatives,
! including OLAM:  

   !----------------------------------------------------------------------------
   ! Copyright (C) 1991-2006  ; All Rights Reserved ; Colorado State University; 
   ! Colorado State University Research Foundation ; ATMET, LLC 

   ! This software is free software; you can redistribute it and/or modify it 
   ! under the terms of the GNU General Public License as published by the Free
   ! Software Foundation; either version 2 of the License, or (at your option)
   ! any later version. 

   ! This software is distributed in the hope that it will be useful, but
   ! WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
   ! or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
   ! for more details.
 
   ! You should have received a copy of the GNU General Public License along
   ! with this program; if not, write to the Free Software Foundation, Inc.,
   ! 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA 
   ! (http://www.gnu.org/licenses/gpl.html) 
   !----------------------------------------------------------------------------

! It is requested that any scientific publications based on application of OLAM
! include the following acknowledgment:  "OLAM was developed at the 
! Edmund T. Pratt Jr. School of Engineering, Duke University."

! For additional information, including published references, please contact
! the software authors, Robert L. Walko (robert.walko@duke.edu)
! or Roni Avissar (avissar@duke.edu).
!===============================================================================
subroutine harr_radinit()  ! Bob's interface subroutine

use mem_harr, only: ng,nb,nsolb,npsb,nuum,prf,alpha,trf,beta  &
                    ,xp,wght,wlenlo,wlenhi,solar0,ralcs,a0,a1,a2,a3  &
                    ,ulim,npartob,npartg,ncog,ncb  &
                    ,ocoef,bcoef,gcoef,solara,solarb

use mem_radiate, only: maxadd_rad, nadd_rad, zmrad
use micro_coms,  only: gnu
use mem_grid,    only: mza, zm
use misc_coms,   only: io6

implicit none

real :: deltaz

!------------------------------------------------------------------------------
! This subroutine initializes some constants for the Harrington radiation
! scheme.  It is called only once, by the master process or by each 
! compute node process, at the time of model initialization.
!------------------------------------------------------------------------------

! Initialize constants for Harrington s/w and l/w radiation computations

call harr_radinit1(ng,nb,nsolb,npsb,nuum,prf,alpha,trf,beta  &
   ,xp,wght,wlenlo,wlenhi,solar0,ralcs,a0,a1,a2,a3  &
   ,ulim,npartob,npartg,ncog,ncb  &
   ,ocoef,bcoef,gcoef,gnu)

! Compute NADD_RAD, the number of radiation levels to be added above the top
! model prognostic level.  (Added levels will be filled elsewhere with data 
! from Mclatchy soundings.)

if (zm(mza-1) > 25000.) then
   
   ! If model top is above 25 km (roughly 30 mb tropics), do not add any levels.

   nadd_rad = 0
   zmrad = zm(mza-1)

else

   ! If grid top < 25 km, add one or more radiation levels.
   ! Set the top level at 30 km (rougly 15 mb tropics).

   ! Estimate a reasonable value for the height increment between added levels.
   ! Make it approximately the increment between the two highest model levels,
   ! but no less than allowed by maxadd_rad, the maximum number of added levels.
   
   zmrad = 30.e3

   deltaz = max( zm(mza-1) - zm(mza-2), (zmrad - zm(mza-1)) / real(maxadd_rad) )

   nadd_rad = nint( (zmrad - zm(mza-1)) / deltaz )
   nadd_rad = max(nadd_rad,1)
   nadd_rad = min(nadd_rad,maxadd_rad)

endif

return
end subroutine harr_radinit

!===============================================================================

subroutine harr_radinit1(ng,nb,ns,npsb,nuum,prf,alpha,trf,beta,xp,  &
    wght,wlenlo,wlenhi,solar,ralcs,a0,a1,a2,a3,ulim,  &
    npartob,npartg,ncog,ncb,ocoef,bcoef,gcoef,gnu)
    
implicit none

!
!     read input parameters for the radiative transfer scheme
!     Original: J Verlinde March 1993
!
!     Rewritten and modified: J.H. 1995/1996
!
integer, parameter :: mb=8,mg=3,mk=7
real, parameter :: top=1800.,tm=1800./296.
!
!     mb   --   maximum number of bands allowed
!     mg   --   maximum active gasses in scheme (different from
!               the max number of active gasses in any given band)
!     mk   --   maximum number of pseudobands allowed for any gas
!     (H2O - gas 1;  CO2 plus other gases - gas 2;  Ozone - gas 3)
!

integer :: ng,ns,nb,npsb(mg,mb),nuum(mb),npartob,npartg,ncog,ncb
real :: prf(mg,mb),alpha(mg,mk,mb),trf(mg,mb),beta(mg,mk,mb),  &
     xp(mg,mk,mb),wght(mg,mk,mb),  &
     wlenlo(mb),wlenhi(mb),solar(mb),ralcs(mb),  &
     a0(mb),a1(mb),a2(mb),a3(mb),solc,ulim(mg,mb)
real :: ocoef(ncog,mb,npartob),bcoef(ncb,mb,npartob),  &
     gcoef(ncog,mb,npartg)
real, intent(in) :: gnu(7)

integer :: npsb1(3,8),nuum1(8)
real :: prf1(3,8),alpha1(3,7,8),trf1(3,8),beta1(3,7,8),  &
     xp1(3,7,8),wght1(3,7,8),  &
     wlenlo1(8),wlenhi1(8),solar1(8),ralcs1(8),  &
     a01(8),a11(8),a21(8),a31(8),  &
     ulim1(3,8)
real :: ocoef1(5,8,13,2),bcoef1(2,8,13,2)  &
    ,gcoef1(5,8,7,2)

integer :: ib,icog,icb,ig,ik,ngb,ip,icat,ignu,ic1,ic2,ic3,i,is
real :: sunavg,temp

integer :: kkat(13)
data kkat /1,2,3,3,3,4,4,4,5,5,5,6,7/

data ((ocoef1(icog,ib, 1,1),icog=1,5),ib=1,8)/  &
   .53502E+00,  .25038E+00, -.28307E-02,  .17203E+00, -.22921E-01,  &
   .87670E+00, -.14653E-02, -.22618E+01,  .12308E+00, -.92743E-03,  &
   .99986E+00, -.49475E+00, -.19301E+01,  .60422E+00, -.21285E+01,  &
   .51950E+00, -.14330E+00, -.68641E-01, -.14036E+01, -.18877E+01,  &
   .51500E+00, -.75507E+00, -.25649E+00,  .16039E+01, -.19427E+01,  &
   .51014E+00, -.78594E+00, -.29105E+01,  .21041E+00, -.34407E-01,  &
   .50916E+00,  .68010E-01, -.12197E-01, -.19787E+01, -.25706E+01,  &
   .50477E+00,  .26895E+00, -.50744E-01,  .39204E-01, -.65051E-02/
data ((ocoef1(icog,ib, 2,1),icog=1,5),ib=1,8)/  &
   .50246E+00,  .27495E+00, -.86761E-02,  .15534E+00, -.11628E-02,  &
   .77184E+00,  .21225E+01, -.46616E+00,  .20555E+00, -.26433E-03,  &
   .99864E+00,  .29610E+01, -.76111E+00,  .24400E+01, -.21687E+01,  &
   .50373E+00, -.56477E+00, -.18834E+00,  .27339E-01, -.95556E-03,  &
   .50422E+00, -.10954E+01, -.32495E+00,  .26231E+01, -.12727E+01,  &
   .50285E+00,  .11969E+00, -.13137E-01,  .29679E+01, -.35296E+00,  &
   .50250E+00,  .96104E+00, -.36542E+00,  .70906E-01, -.63027E-02,  &
   .50217E+00,  .85769E+00, -.15454E+00,  .35989E-01, -.46899E-02/
data ((ocoef1(icog,ib, 3,1),icog=1,5),ib=1,8)/  &
   .81912E+00,  .89934E-01, -.95431E-02,  .77550E-01, -.64635E-01,  &
   .99919E+00,  .21984E+01, -.26508E+01, -.25080E+01, -.27903E+01,  &
   .99995E+00,  .30675E+00, -.24977E+01, -.40199E+00, -.27822E+01,  &
   .66707E+00,  .21413E+00, -.19915E-01, -.14337E+00, -.39737E+00,  &
   .56341E+00,  .21768E+00, -.47805E-01, -.39828E+00, -.15905E+01,  &
   .56926E+00, -.28098E+00, -.12872E+01,  .35146E+00, -.39604E-01,  &
   .55241E+00, -.21384E+00, -.80189E+00,  .18383E+00, -.36890E-01,  &
   .52773E+00,  .28583E+00, -.69522E-01,  .13389E+00, -.10619E-01/
data ((ocoef1(icog,ib, 4,1),icog=1,5),ib=1,8)/  &
   .88678E+00,  .95164E-01, -.18882E-01, -.79279E-01, -.26767E+01,  &
   .99955E+00, -.95345E+00, -.27709E+01,  .95466E+00, -.27609E+01,  &
   .99995E+00,  .30675E+00, -.24977E+01, -.40199E+00, -.27822E+01,  &
   .65337E+00, -.29702E+00, -.86659E+00,  .24877E+00, -.63270E-02,  &
   .57824E+00,  .23396E+00, -.29201E-01, -.97145E+00, -.20002E+01,  &
   .60869E+00, -.54567E+00, -.16368E+01,  .32509E+00, -.19304E-01,  &
   .57423E+00,  .20302E+00, -.20766E-01, -.52963E+00, -.12669E+01,  &
   .62105E+00,  .31371E+00, -.24231E-01, -.83159E+00, -.28613E+01/
data ((ocoef1(icog,ib, 5,1),icog=1,5),ib=1,8)/  &
   .85390E+00,  .24811E-01, -.59938E-01,  .97555E-01, -.13096E-01,  &
   .99945E+00,  .19216E+00, -.22428E+01, -.38142E+00, -.29731E+01,  &
   .99995E+00,  .30675E+00, -.24977E+01, -.40199E+00, -.27822E+01,  &
   .59890E+00, -.18424E+00, -.35464E+00,  .26931E+00, -.61489E-02,  &
   .56129E+00,  .20357E+00, -.27843E-01, -.52470E+00, -.13424E+01,  &
   .57703E+00, -.27255E+00, -.88195E+00,  .33391E+00, -.23252E-01,  &
   .54768E+00, -.29256E+00, -.61564E+00,  .17540E+00, -.19055E-01,  &
   .58753E+00,  .32487E+00, -.27846E-01, -.80520E+00, -.28833E+01/
data ((ocoef1(icog,ib, 6,1),icog=1,5),ib=1,8)/  &
   .60131E+00,  .14296E+01, -.26494E-01,  .19873E+00, -.25000E-03,  &
   .99253E+00,  .10755E+00, -.23677E-01, -.91790E+00, -.15857E+00,  &
   .99997E+00, -.27667E+01, -.85050E+00, -.15914E+01, -.28007E+01,  &
   .50331E+00,  .20589E+01, -.24104E+01,  .25282E+00, -.27455E-02,  &
   .50003E+00, -.18764E+01, -.39166E-01,  .12933E+00, -.45107E-02,  &
   .49990E+00, -.16069E+01, -.52634E-01,  .11701E+00, -.44343E-02,  &
   .49979E+00, -.23946E+01, -.22541E+01,  .91240E-01, -.41169E-02,  &
   .50012E+00,  .79855E+00, -.16299E+01,  .10765E+00, -.41424E-02/
data ((ocoef1(icog,ib, 7,1),icog=1,5),ib=1,8)/  &
   .61674E+00,  .21205E+01, -.29820E-01,  .23706E+00, -.22189E-03,  &
   .99432E+00,  .13598E+01, -.44456E-01,  .12953E+01, -.14623E+01,  &
   .99998E+00,  .21184E+01, -.92367E-01,  .53557E+00, -.35783E+00,  &
   .50773E+00,  .27353E+01, -.35146E-01,  .25352E+00, -.14164E-02,  &
   .49986E+00, -.23090E+01, -.14960E+01,  .13935E+00, -.41030E-02,  &
   .50042E+00,  .22300E+00, -.41742E-02, -.22590E+01, -.22655E+00,  &
   .50006E+00, -.24034E+01, -.19547E+01,  .14781E+00, -.41357E-02,  &
   .50123E+00, -.29943E+01, -.29163E+01,  .18482E+00, -.32265E-02/
data ((ocoef1(icog,ib, 8,1),icog=1,5),ib=1,8)/  &
   .60552E+00,  .19524E+00, -.91893E-02,  .22414E+00, -.22649E-03,  &
   .99345E+00, -.25853E+01, -.11440E+01,  .76242E+00, -.37469E-01,  &
   .99998E+00,  .96513E+00, -.84851E-01,  .29527E+01, -.11622E+01,  &
   .50469E+00,  .21314E+00, -.17656E-02,  .33216E+00, -.13907E-01,  &
   .50001E+00, -.10675E+01, -.52958E-01,  .11642E+00, -.42821E-02,  &
   .50003E+00, -.18789E+01, -.15510E+01,  .16530E+00, -.47971E-02,  &
   .50005E+00,  .11219E+00, -.44481E-02,  .15162E+00, -.11005E+01,  &
   .50059E+00,  .15061E+00, -.40776E-02, -.19571E+01, -.19943E+01/
data ((ocoef1(icog,ib, 9,1),icog=1,5),ib=1,8)/  &
   .59111E+00,  .19384E+00, -.20844E-03,  .15061E+00, -.50790E-02,  &
   .99242E+00,  .11224E+01, -.50262E+00, -.11549E+01, -.24816E+01,  &
   .99997E+00,  .10566E+01, -.12959E+01,  .29739E+01, -.11756E+01,  &
   .50337E+00,  .96222E+00, -.23622E+00,  .25860E+00, -.28177E-02,  &
   .49988E+00,  .10651E+00, -.42909E-02,  .17299E+01, -.28869E+00,  &
   .49997E+00,  .10887E+01, -.17114E+00,  .12077E+00, -.46694E-02,  &
   .49987E+00,  .94228E-01, -.43431E-02,  .54459E+00, -.17733E+00,  &
   .50023E+00,  .10820E+00, -.41047E-02,  .28487E+01, -.26901E+00/
data ((ocoef1(icog,ib,10,1),icog=1,5),ib=1,8)/  &
   .59885E+00,  .12438E+00, -.37180E-02,  .23437E+00, -.17364E-03,  &
   .99431E+00,  .98683E+00, -.50705E+00, -.15639E+01, -.17063E+01,  &
   .99998E+00,  .24721E+01, -.18065E+01,  .23966E+01, -.14163E+01,  &
   .50310E+00,  .79901E-01, -.51442E-03,  .30571E+00, -.34180E-02,  &
   .49967E+00,  .98947E+00, -.20676E+00,  .13659E+00, -.41175E-02,  &
   .50059E+00,  .25137E+00, -.39431E-01,  .22353E+00, -.42388E-02,  &
   .50008E+00,  .15606E+00, -.43760E-02,  .19032E+01, -.30191E+00,  &
   .50149E+00,  .19328E+00, -.33401E-02,  .16638E+01, -.22392E+00/
data ((ocoef1(icog,ib,11,1),icog=1,5),ib=1,8)/  &
   .59114E+00,  .22326E+00, -.18499E-03,  .13978E+00, -.47685E-02,  &
   .99306E+00,  .24021E+01, -.91627E+00,  .42397E-02, -.39950E-03,  &
   .99998E+00, -.27872E+01, -.17821E+01,  .27485E+01, -.11735E+01,  &
   .50579E+00,  .28701E+01, -.36958E+00,  .27189E+00, -.22118E-02,  &
   .50008E+00,  .27605E+01, -.32258E+00,  .11675E+00, -.44468E-02,  &
   .49984E+00,  .29173E+01, -.27598E+00,  .16851E+00, -.48549E-02,  &
   .50011E+00,  .13082E+00, -.50924E-02,  .10954E+01, -.28379E+00,  &
   .50050E+00,  .73313E+00, -.13738E+00,  .15749E+00, -.42053E-02/
data ((ocoef1(icog,ib,12,1),icog=1,5),ib=1,8)/  &
   .50246E+00,  .27495E+00, -.86761E-02,  .15534E+00, -.11628E-02,  &
   .77184E+00,  .21225E+01, -.46616E+00,  .20555E+00, -.26433E-03,  &
   .99864E+00,  .29610E+01, -.76111E+00,  .24400E+01, -.21687E+01,  &
   .50373E+00, -.56477E+00, -.18834E+00,  .27339E-01, -.95556E-03,  &
   .50422E+00, -.10954E+01, -.32495E+00,  .26231E+01, -.12727E+01,  &
   .50285E+00,  .11969E+00, -.13137E-01,  .29679E+01, -.35296E+00,  &
   .50250E+00,  .96104E+00, -.36542E+00,  .70906E-01, -.63027E-02,  &
   .50217E+00,  .85769E+00, -.15454E+00,  .35989E-01, -.46899E-02/
data ((ocoef1(icog,ib,13,1),icog=1,5),ib=1,8)/  &
   .50246E+00,  .27495E+00, -.86761E-02,  .15534E+00, -.11628E-02,  &
   .77184E+00,  .21225E+01, -.46616E+00,  .20555E+00, -.26433E-03,  &
   .99864E+00,  .29610E+01, -.76111E+00,  .24400E+01, -.21687E+01,  &
   .50373E+00, -.56477E+00, -.18834E+00,  .27339E-01, -.95556E-03,  &
   .50422E+00, -.10954E+01, -.32495E+00,  .26231E+01, -.12727E+01,  &
   .50285E+00,  .11969E+00, -.13137E-01,  .29679E+01, -.35296E+00,  &
   .50250E+00,  .96104E+00, -.36542E+00,  .70906E-01, -.63027E-02,  &
   .50217E+00,  .85769E+00, -.15454E+00,  .35989E-01, -.46899E-02/

data ((ocoef1(icog,ib, 1,2),icog=1,5),ib=1,8)/  &
   .52438E+00,  .24743E+00, -.34625E-02,  .18970E+00, -.29927E-01,  &
   .86366E+00,  .28636E-01, -.22157E+01,  .13490E+00, -.10851E-02,  &
   .99979E+00, -.29198E+01, -.28868E+01,  .18287E+01, -.24188E+01,  &
   .49952E+00,  .31959E-01, -.12977E-02, -.39239E+00, -.15836E+00,  &
   .50631E+00, -.43099E+00, -.21672E+00,  .29032E-01, -.42007E-02,  &
   .50446E+00,  .18046E+00, -.64406E-01,  .25188E-01, -.51384E-02,  &
   .50732E+00, -.20182E+01, -.27063E+01,  .61180E-01, -.13667E-01,  &
   .50425E+00,  .37006E-01, -.90129E-02,  .29626E+00, -.79197E-01/
data ((ocoef1(icog,ib, 2,2),icog=1,5),ib=1,8)/  &
   .50131E+00,  .27610E+00, -.10296E-01,  .13972E+00, -.14901E-02,  &
   .75380E+00,  .21442E+00, -.29393E-03,  .28470E+00, -.20896E+00,  &
   .99820E+00, -.29567E+01, -.74150E+00,  .13607E+01, -.59503E+00,  &
   .50496E+00, -.29518E+01, -.87080E+00, -.21751E+01, -.40022E+00,  &
   .50338E+00, -.29793E+01, -.93078E+00, -.24036E+01, -.51415E+00,  &
   .50201E+00,  .27620E+01, -.52448E+00,  .16007E+00, -.18536E-01,  &
   .50190E+00,  .23171E+01, -.36209E+00,  .19042E-01, -.18868E-02,  &
   .50138E+00,  .18239E+00, -.21161E-01,  .21794E+01, -.45835E+00/
data ((ocoef1(icog,ib, 3,2),icog=1,5),ib=1,8)/  &
   .81133E+00,  .69349E-01, -.10299E+00,  .10795E+00, -.14262E-01,  &
   .99895E+00, -.20536E+01, -.25147E+01,  .17417E+01, -.23481E+01,  &
   .99995E+00, -.40557E+00, -.26487E+01,  .36965E+00, -.25690E+01,  &
   .64454E+00,  .23063E+00, -.22302E-01, -.15667E+00, -.67213E+00,  &
   .55383E+00, -.70583E+00, -.26693E+01,  .22107E+00, -.59426E-01,  &
   .61802E+00, -.34675E-01,  .61173E-02,  .33720E+00, -.60327E-01,  &
   .54433E+00, -.19635E+00, -.91939E+00,  .19782E+00, -.48602E-01,  &
   .52284E+00,  .11900E+00, -.13524E-01,  .31394E+00, -.91911E-01/
data ((ocoef1(icog,ib, 4,2),icog=1,5),ib=1,8)/  &
   .87233E+00,  .99292E-01, -.19262E-01,  .17117E-01, -.19705E+00,  &
   .99953E+00, -.16322E+00, -.29089E+01,  .84448E-01, -.21449E+01,  &
   .99995E+00,  .30675E+00, -.24977E+01, -.40199E+00, -.27822E+01,  &
   .63934E+00,  .26529E+00, -.85349E-02, -.21627E+00, -.78589E+00,  &
   .56400E+00,  .24216E+00, -.36364E-01, -.14425E+01, -.29162E+01,  &
   .58293E+00,  .34643E+00, -.23927E-01, -.13853E+01, -.29851E+01,  &
   .55724E+00,  .21398E+00, -.24772E-01, -.80610E+00, -.20961E+01,  &
   .54695E+00,  .21108E+00, -.11528E-01,  .18701E+00, -.59256E-01/
data ((ocoef1(icog,ib, 5,2),icog=1,5),ib=1,8)/  &
   .84401E+00,  .23811E-01, -.14551E+00,  .10861E+00, -.15419E-01,  &
   .99932E+00, -.69472E+00, -.23245E+01,  .67102E+00, -.22771E+01,  &
   .99995E+00,  .30675E+00, -.24977E+01, -.40199E+00, -.27822E+01,  &
   .99963E+00, -.16618E+00,  .53578E-02, -.22109E+00, -.10289E+01,  &
   .55164E+00,  .21466E+00, -.36566E-01, -.60797E+00, -.18416E+01,  &
   .56015E+00, -.32017E+00, -.13440E+01,  .34671E+00, -.29540E-01,  &
   .54148E+00, -.25292E+00, -.72859E+00,  .18398E+00, -.25635E-01,  &
   .53004E+00,  .22726E+00, -.60634E-01,  .16822E+00, -.12088E-01/
data ((ocoef1(icog,ib, 6,2),icog=1,5),ib=1,8)/  &
   .56528E+00,  .20898E+00, -.23015E-03,  .18113E+00, -.66666E-02,  &
   .99010E+00,  .51630E-01, -.10120E-01, -.23299E+01, -.68731E-01,  &
   .99997E+00,  .16679E+01, -.88622E-01,  .18730E+01, -.63071E+00,  &
   .50148E+00,  .13349E+00, -.25006E+01,  .23413E+00, -.32827E-02,  &
   .49998E+00, -.28983E+01, -.40580E-01,  .12123E+00, -.52029E-02,  &
   .49984E+00,  .79700E-01, -.42625E-02, -.63772E+00, -.92501E+00,  &
   .49990E+00,  .73620E-01, -.41556E-02,  .22763E+01, -.24832E+01,  &
   .49997E+00,  .26563E+01, -.15949E+01,  .78663E-01, -.42985E-02/
data ((ocoef1(icog,ib, 7,2),icog=1,5),ib=1,8)/  &
   .56310E+00,  .25170E+00, -.17797E-03,  .88188E-01, -.24857E-02,  &
   .99255E+00,  .26672E+01, -.47333E-01, -.13995E+01, -.24663E+01,  &
   .99998E+00,  .18766E+01, -.18795E+01, -.85618E-01, -.19384E+01,  &
   .50248E+00,  .89183E-01, -.93480E-03,  .26905E+00, -.41338E-02,  &
   .50020E+00, -.11036E+01, -.30119E-01,  .18127E+00, -.59058E-02,  &
   .50014E+00, -.19260E+01, -.13113E+01,  .17751E+00, -.49112E-02,  &
   .50003E+00,  .11324E+00, -.45075E-02, -.17264E+00, -.12135E+01,  &
   .50038E+00,  .15476E+00, -.39471E-02,  .14488E+01, -.13307E+01/
data ((ocoef1(icog,ib, 8,2),icog=1,5),ib=1,8)/  &
   .58277E+00,  .24129E+00, -.25302E-03,  .29461E+01, -.32957E-01,  &
   .99189E+00,  .60466E+00, -.35096E-01, -.26453E+01, -.83121E+00,  &
   .99997E+00, -.27667E+01, -.85050E+00, -.15914E+01, -.28007E+01,  &
   .50333E+00,  .22535E+00, -.24322E-02,  .29507E+01, -.37696E-01,  &
   .49987E+00, -.19453E+01, -.29072E+01,  .89687E-01, -.42578E-02,  &
   .50009E+00, -.25865E+00, -.20818E+01,  .12573E+00, -.53493E-02,  &
   .49994E+00, -.99434E+00, -.13396E+01,  .85224E-01, -.46679E-02,  &
   .50013E+00,  .11770E+00, -.45465E-02, -.85177E+00, -.27648E+01/
data ((ocoef1(icog,ib, 9,2),icog=1,5),ib=1,8)/  &
   .56278E+00,  .20772E+00, -.21767E-03,  .15660E+00, -.52773E-02,  &
   .99044E+00, -.29442E+01, -.13689E+01,  .24217E+01, -.55470E+00,  &
   .99997E+00,  .96718E+00, -.25363E+01,  .24388E+01, -.11205E+01,  &
   .50139E+00,  .23373E+00, -.32921E-02,  .15407E+01, -.26574E+00,  &
   .49978E+00,  .84373E-01, -.38509E-02,  .62946E+00, -.19821E+00,  &
   .49997E+00,  .14601E+01, -.20115E+00,  .83662E-01, -.43680E-02,  &
   .49976E+00,  .25618E+01, -.32259E+00,  .71141E-01, -.41767E-02,  &
   .49958E+00,  .73392E-01, -.39063E-02,  .23186E+01, -.25438E+00/
data ((ocoef1(icog,ib,10,2),icog=1,5),ib=1,8)/  &
   .57280E+00,  .13298E+00, -.53572E-02,  .25483E+00, -.20239E-03,  &
   .99265E+00,  .11151E+00, -.22915E+01,  .29903E+01, -.60607E+00,  &
   .99998E+00, -.17678E+01, -.15604E+01,  .95765E+00, -.10485E+01,  &
   .50455E+00,  .19910E+01, -.28865E+00,  .27178E+00, -.20250E-02,  &
   .49994E+00,  .24198E+01, -.29927E+00,  .11348E+00, -.44506E-02,  &
   .50015E+00,  .27335E+00, -.33081E-01,  .17208E+00, -.49023E-02,  &
   .50019E+00,  .11259E+00, -.44642E-02,  .24168E+00, -.70501E-01,  &
   .50052E+00,  .15972E+00, -.40628E-02,  .47993E+00, -.93927E-01/
data ((ocoef1(icog,ib,11,2),icog=1,5),ib=1,8)/  &
   .56714E+00,  .14168E+00, -.60114E-02,  .23864E+00, -.20326E-03,  &
   .99203E+00,  .27487E+00, -.35208E+00, -.11985E+00, -.22339E+01,  &
   .99997E+00,  .28339E+01, -.12505E+01, -.28704E+01, -.29825E+01,  &
   .50298E+00,  .60429E+00, -.18340E+00,  .25075E+00, -.25888E-02,  &
   .50004E+00,  .86340E-01, -.44164E-02,  .22731E+01, -.29232E+00,  &
   .49986E+00,  .12787E+00, -.54055E-02,  .24263E+01, -.24682E+00,  &
   .49980E+00,  .85550E-01, -.43925E-02,  .10047E+01, -.22929E+00,  &
   .50024E+00,  .10162E+01, -.16810E+00,  .11731E+00, -.45331E-02/
data ((ocoef1(icog,ib,12,2),icog=1,5),ib=1,8)/  &
   .50131E+00,  .27610E+00, -.10296E-01,  .13972E+00, -.14901E-02,  &
   .75380E+00,  .21442E+00, -.29393E-03,  .28470E+00, -.20896E+00,  &
   .99820E+00, -.29567E+01, -.74150E+00,  .13607E+01, -.59503E+00,  &
   .50496E+00, -.29518E+01, -.87080E+00, -.21751E+01, -.40022E+00,  &
   .50338E+00, -.29793E+01, -.93078E+00, -.24036E+01, -.51415E+00,  &
   .50201E+00,  .27620E+01, -.52448E+00,  .16007E+00, -.18536E-01,  &
   .50190E+00,  .23171E+01, -.36209E+00,  .19042E-01, -.18868E-02,  &
   .50138E+00,  .18239E+00, -.21161E-01,  .21794E+01, -.45835E+00/
data ((ocoef1(icog,ib,13,2),icog=1,5),ib=1,8)/  &
   .50131E+00,  .27610E+00, -.10296E-01,  .13972E+00, -.14901E-02,  &
   .75380E+00,  .21442E+00, -.29393E-03,  .28470E+00, -.20896E+00,  &
   .99820E+00, -.29567E+01, -.74150E+00,  .13607E+01, -.59503E+00,  &
   .50496E+00, -.29518E+01, -.87080E+00, -.21751E+01, -.40022E+00,  &
   .50338E+00, -.29793E+01, -.93078E+00, -.24036E+01, -.51415E+00,  &
   .50201E+00,  .27620E+01, -.52448E+00,  .16007E+00, -.18536E-01,  &
   .50190E+00,  .23171E+01, -.36209E+00,  .19042E-01, -.18868E-02,  &
   .50138E+00,  .18239E+00, -.21161E-01,  .21794E+01, -.45835E+00/

data ((bcoef1(icb,ib, 1,1),icb=1,2),ib=1,8)/  &
   .34438E-11,  .19866E+01,  &
   .34388E-11,  .19863E+01,  &
   .33441E-11,  .19905E+01,  &
   .32585E-11,  .20009E+01,  &
   .31041E-11,  .20064E+01,  &
   .30466E-11,  .20084E+01,  &
   .29263E-11,  .20152E+01,  &
   .31930E-11,  .20001E+01/
data ((bcoef1(icb,ib, 2,1),icb=1,2),ib=1,8)/  &
   .32688E-11,  .19956E+01,  &
   .32122E-11,  .19975E+01,  &
   .31842E-11,  .19985E+01,  &
   .35364E-11,  .19873E+01,  &
   .33857E-11,  .19920E+01,  &
   .34428E-11,  .19898E+01,  &
   .33757E-11,  .19922E+01,  &
   .33885E-11,  .19916E+01/
data ((bcoef1(icb,ib, 3,1),icb=1,2),ib=1,8)/  &
   .58083E-11,  .17438E+01,  &
   .55947E-11,  .17484E+01,  &
   .52795E-11,  .17598E+01,  &
   .34851E-11,  .19093E+01,  &
   .54453E-11,  .17832E+01,  &
   .51652E-11,  .17924E+01,  &
   .41981E-11,  .18363E+01,  &
   .58427E-11,  .17551E+01/
data ((bcoef1(icb,ib, 4,1),icb=1,2),ib=1,8)/  &
   .18179E-10,  .14989E+01,  &
   .17374E-10,  .15048E+01,  &
   .16215E-10,  .15185E+01,  &
   .89107E-11,  .17122E+01,  &
   .16339E-10,  .15501E+01,  &
   .15479E-10,  .15601E+01,  &
   .11922E-10,  .16147E+01,  &
   .18206E-10,  .15130E+01/
data ((bcoef1(icb,ib, 5,1),icb=1,2),ib=1,8)/  &
   .48019E-11,  .17647E+01,  &
   .45784E-11,  .17712E+01,  &
   .42466E-11,  .17864E+01,  &
   .20760E-11,  .20058E+01,  &
   .41150E-11,  .18264E+01,  &
   .38476E-11,  .18393E+01,  &
   .29229E-11,  .18972E+01,  &
   .47067E-11,  .17835E+01/
data ((bcoef1(icb,ib, 6,1),icb=1,2),ib=1,8)/  &
   .17941E-10,  .14738E+01,  &
   .17785E-10,  .14748E+01,  &
   .17697E-10,  .14754E+01,  &
   .20203E-10,  .14601E+01,  &
   .18985E-10,  .14673E+01,  &
   .18692E-10,  .14691E+01,  &
   .18672E-10,  .14692E+01,  &
   .18378E-10,  .14710E+01/
data ((bcoef1(icb,ib, 7,1),icb=1,2),ib=1,8)/  &
   .10442E-10,  .16066E+01,  &
   .10358E-10,  .16075E+01,  &
   .10310E-10,  .16081E+01,  &
   .11679E-10,  .15936E+01,  &
   .11018E-10,  .16004E+01,  &
   .10855E-10,  .16021E+01,  &
   .10848E-10,  .16022E+01,  &
   .10687E-10,  .16039E+01/
data ((bcoef1(icb,ib, 8,1),icb=1,2),ib=1,8)/  &
   .16713E-11,  .19739E+01,  &
   .16595E-11,  .19747E+01,  &
   .16533E-11,  .19752E+01,  &
   .18340E-11,  .19632E+01,  &
   .17481E-11,  .19687E+01,  &
   .17252E-11,  .19703E+01,  &
   .17254E-11,  .19702E+01,  &
   .17030E-11,  .19717E+01/
data ((bcoef1(icb,ib, 9,1),icb=1,2),ib=1,8)/  &
   .15886E-10,  .14881E+01,  &
   .15600E-10,  .14903E+01,  &
   .15437E-10,  .14915E+01,  &
   .17739E-10,  .14753E+01,  &
   .17223E-10,  .14787E+01,  &
   .17031E-10,  .14800E+01,  &
   .16499E-10,  .14837E+01,  &
   .16659E-10,  .14826E+01/
data ((bcoef1(icb,ib,10,1),icb=1,2),ib=1,8)/  &
   .11227E-10,  .15981E+01,  &
   .10998E-10,  .16005E+01,  &
   .10864E-10,  .16019E+01,  &
   .12060E-10,  .15899E+01,  &
   .12088E-10,  .15895E+01,  &
   .11974E-10,  .15906E+01,  &
   .11453E-10,  .15958E+01,  &
   .11794E-10,  .15924E+01/
data ((bcoef1(icb,ib,11,1),icb=1,2),ib=1,8)/  &
   .18718E-11,  .19607E+01,  &
   .18382E-11,  .19628E+01,  &
   .18167E-11,  .19642E+01,  &
   .19564E-11,  .19556E+01,  &
   .19861E-11,  .19538E+01,  &
   .19683E-11,  .19549E+01,  &
   .18873E-11,  .19598E+01,  &
   .19525E-11,  .19558E+01/
data ((bcoef1(icb,ib,12,1),icb=1,2),ib=1,8)/  &
   .32688E-11,  .19956E+01,  &
   .32122E-11,  .19975E+01,  &
   .31842E-11,  .19985E+01,  &
   .35364E-11,  .19873E+01,  &
   .33857E-11,  .19920E+01,  &
   .34428E-11,  .19898E+01,  &
   .33757E-11,  .19922E+01,  &
   .33885E-11,  .19916E+01/
data ((bcoef1(icb,ib,13,1),icb=1,2),ib=1,8)/  &
   .32688E-11,  .19956E+01,  &
   .32122E-11,  .19975E+01,  &
   .31842E-11,  .19985E+01,  &
   .35364E-11,  .19873E+01,  &
   .33857E-11,  .19920E+01,  &
   .34428E-11,  .19898E+01,  &
   .33757E-11,  .19922E+01,  &
   .33885E-11,  .19916E+01/

data ((bcoef1(icb,ib, 1,2),icb=1,2),ib=1,8)/  &
   .10388E-10,  .19854E+01,  &
   .10167E-10,  .19885E+01,  &
   .98859E-11,  .19927E+01,  &
   .99564E-11,  .19966E+01,  &
   .95604E-11,  .20013E+01,  &
   .92513E-11,  .20056E+01,  &
   .89816E-11,  .20107E+01,  &
   .97704E-11,  .19962E+01/
data ((bcoef1(icb,ib, 2,2),icb=1,2),ib=1,8)/  &
   .96951E-11,  .19969E+01,  &
   .95759E-11,  .19982E+01,  &
   .95177E-11,  .19989E+01,  &
   .10482E-10,  .19885E+01,  &
   .10072E-10,  .19928E+01,  &
   .10166E-10,  .19916E+01,  &
   .10082E-10,  .19926E+01,  &
   .99843E-11,  .19936E+01/
data ((bcoef1(icb,ib, 3,2),icb=1,2),ib=1,8)/  &
   .20024E-10,  .16731E+01,  &
   .18971E-10,  .16826E+01,  &
   .18042E-10,  .16927E+01,  &
   .15730E-10,  .17736E+01,  &
   .20721E-10,  .16862E+01,  &
   .20019E-10,  .16905E+01,  &
   .16809E-10,  .17295E+01,  &
   .21236E-10,  .16694E+01/
data ((bcoef1(icb,ib, 4,2),icb=1,2),ib=1,8)/  &
   .44143E-10,  .15092E+01,  &
   .41297E-10,  .15213E+01,  &
   .38861E-10,  .15337E+01,  &
   .29552E-10,  .16484E+01,  &
   .44903E-10,  .15278E+01,  &
   .43528E-10,  .15316E+01,  &
   .34547E-10,  .15827E+01,  &
   .47297E-10,  .15043E+01/
data ((bcoef1(icb,ib, 5,2),icb=1,2),ib=1,8)/  &
   .12295E-10,  .17893E+01,  &
   .11473E-10,  .18021E+01,  &
   .10745E-10,  .18156E+01,  &
   .75080E-11,  .19480E+01,  &
   .12090E-10,  .18150E+01,  &
   .11598E-10,  .18212E+01,  &
   .90997E-11,  .18752E+01,  &
   .12992E-10,  .17872E+01/
data ((bcoef1(icb,ib, 6,2),icb=1,2),ib=1,8)/  &
   .42978E-10,  .14781E+01,  &
   .42685E-10,  .14789E+01,  &
   .42514E-10,  .14793E+01,  &
   .47449E-10,  .14666E+01,  &
   .45052E-10,  .14726E+01,  &
   .44365E-10,  .14744E+01,  &
   .44464E-10,  .14741E+01,  &
   .43801E-10,  .14759E+01/
data ((bcoef1(icb,ib, 7,2),icb=1,2),ib=1,8)/  &
   .26517E-10,  .16100E+01,  &
   .26350E-10,  .16107E+01,  &
   .26252E-10,  .16112E+01,  &
   .28994E-10,  .15997E+01,  &
   .27688E-10,  .16050E+01,  &
   .27286E-10,  .16067E+01,  &
   .27348E-10,  .16064E+01,  &
   .26969E-10,  .16080E+01/
data ((bcoef1(icb,ib, 8,2),icb=1,2),ib=1,8)/  &
   .47812E-11,  .19789E+01,  &
   .47558E-11,  .19795E+01,  &
   .47427E-11,  .19798E+01,  &
   .51322E-11,  .19707E+01,  &
   .49508E-11,  .19749E+01,  &
   .48910E-11,  .19763E+01,  &
   .49022E-11,  .19760E+01,  &
   .48465E-11,  .19773E+01/
data ((bcoef1(icb,ib, 9,2),icb=1,2),ib=1,8)/  &
   .39775E-10,  .14872E+01,  &
   .39221E-10,  .14888E+01,  &
   .38924E-10,  .14897E+01,  &
   .45126E-10,  .14725E+01,  &
   .42953E-10,  .14782E+01,  &
   .42418E-10,  .14797E+01,  &
   .41699E-10,  .14817E+01,  &
   .41391E-10,  .14825E+01/
data ((bcoef1(icb,ib,10,2),icb=1,2),ib=1,8)/  &
   .28126E-10,  .16031E+01,  &
   .27668E-10,  .16050E+01,  &
   .27422E-10,  .16061E+01,  &
   .30896E-10,  .15922E+01,  &
   .30342E-10,  .15943E+01,  &
   .30081E-10,  .15953E+01,  &
   .29174E-10,  .15989E+01,  &
   .29387E-10,  .15980E+01/
data ((bcoef1(icb,ib,11,2),icb=1,2),ib=1,8)/  &
   .52306E-11,  .19684E+01,  &
   .51549E-11,  .19701E+01,  &
   .51143E-11,  .19710E+01,  &
   .55947E-11,  .19606E+01,  &
   .55687E-11,  .19611E+01,  &
   .55284E-11,  .19620E+01,  &
   .53706E-11,  .19654E+01,  &
   .54315E-11,  .19640E+01/
data ((bcoef1(icb,ib,12,2),icb=1,2),ib=1,8)/  &
   .96951E-11,  .19969E+01,  &
   .95759E-11,  .19982E+01,  &
   .95177E-11,  .19989E+01,  &
   .10482E-10,  .19885E+01,  &
   .10072E-10,  .19928E+01,  &
   .10166E-10,  .19916E+01,  &
   .10082E-10,  .19926E+01,  &
   .99843E-11,  .19936E+01/
data ((bcoef1(icb,ib,13,2),icb=1,2),ib=1,8)/  &
   .96951E-11,  .19969E+01,  &
   .95759E-11,  .19982E+01,  &
   .95177E-11,  .19989E+01,  &
   .10482E-10,  .19885E+01,  &
   .10072E-10,  .19928E+01,  &
   .10166E-10,  .19916E+01,  &
   .10082E-10,  .19926E+01,  &
   .99843E-11,  .19936E+01/

data ((gcoef1(icog,ib, 1,1),icog=1,5),ib=1,8)/  &
   .85004E+00, -.23070E+01, -.29873E+01, -.72250E-01, -.13559E-01,  &
   .85583E+00, -.39371E-01, -.53225E-01, -.10528E+01, -.27121E+01,  &
   .89734E+00, -.36884E+00, -.28945E+00,  .29238E+00, -.46323E+00,  &
   .86046E+00, -.23999E+01, -.16198E+01, -.37815E+00, -.69269E-01,  &
   .91260E+00, -.19434E+00, -.77833E-01, -.27120E+01, -.14493E+01,  &
   .97134E+00, -.13436E+01, -.85240E+00, -.15603E+00, -.60289E-01,  &
   .96978E+00, -.27259E+01, -.14606E+01, -.13315E+00, -.55870E-01,  &
   .92327E+00, -.62450E+00, -.51258E+00, -.15262E+00, -.57748E-01/
data ((gcoef1(icog,ib, 2,1),icog=1,5),ib=1,8)/  &
   .85256E+00, -.28551E+01, -.11673E+01, -.87298E-01, -.12332E-01,  &
   .85579E+00, -.58954E+00, -.34927E+00, -.29288E+01, -.58065E+00,  &
   .89934E+00, -.25957E+01, -.51369E+00,  .16269E+01, -.95447E+00,  &
   .86130E+00, -.28143E+00, -.32093E-01,  .15132E+01, -.10311E+01,  &
   .91312E+00, -.11166E+01, -.24026E+00,  .24531E+01, -.64002E+00,  &
   .97157E+00, -.23646E+01, -.32910E+00,  .20974E+01, -.29852E+01,  &
   .97053E+00, -.12651E+00, -.44707E-01, -.27114E+01, -.18854E+01,  &
   .92363E+00, -.93127E+00, -.23277E+00,  .12683E+01, -.10108E+01/
data ((gcoef1(icog,ib, 3,1),icog=1,5),ib=1,8)/  &
   .84652E+00, -.14792E+00, -.49920E+00, -.82844E-01, -.32344E-01,  &
   .86490E+00, -.24280E-01, -.27182E-01, -.10035E+00, -.25304E+00,  &
   .90438E+00, -.57383E-01, -.15491E+00, -.10861E+00, -.94753E+00,  &
   .88805E+00, -.81787E+00, -.17892E+00, -.17723E+00, -.23314E-01,  &
   .89817E+00, -.18697E+00, -.74539E-01, -.94747E+00, -.42580E+00,  &
   .96632E+00, -.17709E+00, -.72026E-01, -.11097E+01, -.70518E+00,  &
   .96598E+00, -.15428E+00, -.70893E-01, -.11102E+01, -.59795E+00,  &
   .91847E+00, -.19007E+00, -.88016E-01, -.10278E+01, -.11070E+01/
data ((gcoef1(icog,ib, 4,1),icog=1,5),ib=1,8)/  &
   .86492E+00,  .14010E+01, -.18127E+01, -.12305E+01, -.34913E-01,  &
   .86413E+00, -.89078E-02, -.40056E+00, -.11315E+01, -.66225E+00,  &
   .90703E+00, -.23060E+01, -.14108E+00,  .16011E+01, -.10697E+01,  &
   .88834E+00, -.10417E+01, -.39147E+00, -.52226E+00, -.29662E-01,  &
   .90180E+00, -.22150E+01, -.26668E+01, -.20787E+01, -.52022E-01,  &
   .96995E+00, -.40797E+00, -.63701E-01, -.48035E+00, -.45402E-01,  &
   .97021E+00, -.22262E+00, -.38171E-01, -.27977E+01, -.75078E+00,  &
   .92271E+00, -.27473E+01, -.60188E-01, -.38875E+00, -.18987E+01/
data ((gcoef1(icog,ib, 5,1),icog=1,5),ib=1,8)/  &
   .86542E+00, -.10243E+00, -.14467E-01,  .18090E+01, -.88272E+00,  &
   .86609E+00, -.33357E+00, -.23570E+01, -.14838E+01, -.40147E+00,  &
   .90692E+00, -.21077E+01, -.67812E+00, -.20880E+01, -.51111E+00,  &
   .88844E+00,  .24039E+01, -.14473E+00, -.12302E+01, -.39473E-01,  &
   .90167E+00, -.18168E+01, -.29193E+00, -.79883E+00, -.43798E+00,  &
   .96988E+00, -.30375E+00, -.12232E+00,  .16400E+01, -.24678E+01,  &
   .97021E+00,  .17677E+01, -.66322E+00, -.16524E+01, -.29933E+00,  &
   .92259E+00, -.16496E+01, -.39920E+00, -.68006E+00, -.26008E+00/
data ((gcoef1(icog,ib, 6,1),icog=1,5),ib=1,8)/  &
   .86542E+00, -.10243E+00, -.14467E-01,  .18090E+01, -.88272E+00,  &
   .86609E+00, -.33357E+00, -.23570E+01, -.14838E+01, -.40147E+00,  &
   .90692E+00, -.21077E+01, -.67812E+00, -.20880E+01, -.51111E+00,  &
   .88844E+00,  .24039E+01, -.14473E+00, -.12302E+01, -.39473E-01,  &
   .90167E+00, -.18168E+01, -.29193E+00, -.79883E+00, -.43798E+00,  &
   .96988E+00, -.30375E+00, -.12232E+00,  .16400E+01, -.24678E+01,  &
   .97021E+00,  .17677E+01, -.66322E+00, -.16524E+01, -.29933E+00,  &
   .92259E+00, -.16496E+01, -.39920E+00, -.68006E+00, -.26008E+00/
data ((gcoef1(icog,ib, 7,1),icog=1,5),ib=1,8)/  &
   .86542E+00, -.10243E+00, -.14467E-01,  .18090E+01, -.88272E+00,  &
   .86609E+00, -.33357E+00, -.23570E+01, -.14838E+01, -.40147E+00,  &
   .90692E+00, -.21077E+01, -.67812E+00, -.20880E+01, -.51111E+00,  &
   .88844E+00,  .24039E+01, -.14473E+00, -.12302E+01, -.39473E-01,  &
   .90167E+00, -.18168E+01, -.29193E+00, -.79883E+00, -.43798E+00,  &
   .96988E+00, -.30375E+00, -.12232E+00,  .16400E+01, -.24678E+01,  &
   .97021E+00,  .17677E+01, -.66322E+00, -.16524E+01, -.29933E+00,  &
   .92259E+00, -.16496E+01, -.39920E+00, -.68006E+00, -.26008E+00/

data ((gcoef1(icog,ib, 1,2),icog=1,5),ib=1,8)/  &
   .85428E+00, -.22824E+00, -.32737E+00, -.57281E-01, -.11823E-01,  &
   .85677E+00, -.15857E+01, -.54104E+00,  .15231E+01, -.63857E+00,  &
   .89818E+00, -.25075E+00, -.24189E+01, -.36307E-01, -.12669E+00,  &
   .85962E+00, -.46391E+00, -.11261E+01, -.62972E+00, -.15502E+00,  &
   .91238E+00, -.15845E+01, -.18563E+01, -.26871E+00, -.14749E+00,  &
   .97117E+00, -.15075E+01, -.15209E+01, -.15929E+00, -.87209E-01,  &
   .96977E+00, -.15247E+00, -.91890E-01, -.29261E+01, -.27680E+01,  &
   .92298E+00, -.16593E+00, -.88094E-01, -.81450E+00, -.89494E+00/
data ((gcoef1(icog,ib, 2,2),icog=1,5),ib=1,8)/  &
   .85485E+00,  .16591E+01, -.13538E+01, -.10826E+00, -.20853E-01,  &
   .85877E+00, -.40050E+00, -.11467E+01, -.22559E+01, -.82045E+00,  &
   .89872E+00, -.48360E+00, -.88272E+00, -.27211E+01, -.10538E+01,  &
   .86142E+00, -.40265E+00, -.68366E-01, -.90545E+00, -.82210E+00,  &
   .91314E+00, -.28791E+00, -.96056E-01,  .16624E+01, -.78796E+00,  &
   .97158E+00, -.17553E+01, -.56189E+00, -.11158E+00, -.14263E+01,  &
   .97046E+00, -.15502E+00, -.10433E+00, -.14530E+01, -.95320E+00,  &
   .92357E+00, -.50976E-01, -.13639E+01, -.21887E+00, -.13027E+00/
data ((gcoef1(icog,ib, 3,2),icog=1,5),ib=1,8)/  &
   .83645E+00, -.10657E+00, -.88227E-01, -.57680E+00, -.29838E+01,  &
   .86231E+00, -.94684E-01, -.28902E+00, -.17472E-01, -.48864E-01,  &
   .90517E+00, -.12820E+00, -.97785E+00, -.30679E-01, -.13764E+00,  &
  -.10789E+00, -.88173E+00, -.17572E+00,  .92677E+00,  .76584E-03,  &
   .90332E+00, -.98943E+00, -.45205E+00, -.75741E-01, -.51938E-01,  &
   .96769E+00, -.14666E+00, -.86753E-01, -.11482E+01, -.95418E+00,  &
   .96824E+00, -.11279E+01, -.76324E+00, -.11148E+00, -.78466E-01,  &
   .91817E+00, -.18217E+00, -.12214E+00, -.11635E+01, -.16800E+01/
data ((gcoef1(icog,ib, 4,2),icog=1,5),ib=1,8)/  &
   .86761E+00, -.25078E+00, -.33773E-01, -.53825E+00, -.92074E+00,  &
   .86637E+00,  .22550E+01, -.19380E+00, -.13276E+01, -.93580E-01,  &
   .90660E+00,  .15513E+01, -.13785E+01, -.18898E+00, -.76385E+00,  &
   .88836E+00, -.24290E+01, -.76275E-01,  .14582E+01, -.16409E+01,  &
   .90184E+00, -.16807E+00, -.64505E-01, -.20434E+01, -.12036E+01,  &
   .96993E+00,  .15686E+01, -.15885E+01, -.87090E+00, -.91118E-01,  &
   .97029E+00, -.19457E+01, -.10327E+00,  .12140E+01, -.16934E+01,  &
   .92260E+00, -.20223E+01, -.12144E+01, -.13366E+01, -.10022E+00/
data ((gcoef1(icog,ib, 5,2),icog=1,5),ib=1,8)/  &
   .86750E+00, -.11043E+00, -.22250E-01, -.26821E+01, -.16937E+01,  &
   .86628E+00,  .40213E+00, -.13349E+01, -.29534E+01, -.86907E+00,  &
   .90633E+00,  .10585E+01, -.20870E+01, -.11491E+01, -.80802E+00,  &
   .88830E+00, -.12621E+00, -.29155E-01, -.46961E+00, -.13071E+00,  &
   .90172E+00, -.71294E+00, -.45464E+00, -.11305E+00, -.65795E-01,  &
   .96998E+00,  .18700E+01, -.22688E+01, -.15064E+00, -.63017E-01,  &
   .97035E+00, -.17586E+00, -.10184E+00,  .62649E+00, -.18503E+01,  &
   .92260E+00,  .10520E+01, -.14577E+01, -.19297E+00, -.11617E+00/
data ((gcoef1(icog,ib, 6,2),icog=1,5),ib=1,8)/  &
   .86750E+00, -.11043E+00, -.22250E-01, -.26821E+01, -.16937E+01,  &
   .86628E+00,  .40213E+00, -.13349E+01, -.29534E+01, -.86907E+00,  &
   .90633E+00,  .10585E+01, -.20870E+01, -.11491E+01, -.80802E+00,  &
   .88830E+00, -.12621E+00, -.29155E-01, -.46961E+00, -.13071E+00,  &
   .90172E+00, -.71294E+00, -.45464E+00, -.11305E+00, -.65795E-01,  &
   .96998E+00,  .18700E+01, -.22688E+01, -.15064E+00, -.63017E-01,  &
   .97035E+00, -.17586E+00, -.10184E+00,  .62649E+00, -.18503E+01,  &
   .92260E+00,  .10520E+01, -.14577E+01, -.19297E+00, -.11617E+00/
data ((gcoef1(icog,ib, 7,2),icog=1,5),ib=1,8)/  &
   .86750E+00, -.11043E+00, -.22250E-01, -.26821E+01, -.16937E+01,  &
   .86628E+00,  .40213E+00, -.13349E+01, -.29534E+01, -.86907E+00,  &
   .90633E+00,  .10585E+01, -.20870E+01, -.11491E+01, -.80802E+00,  &
   .88830E+00, -.12621E+00, -.29155E-01, -.46961E+00, -.13071E+00,  &
   .90172E+00, -.71294E+00, -.45464E+00, -.11305E+00, -.65795E-01,  &
   .96998E+00,  .18700E+01, -.22688E+01, -.15064E+00, -.63017E-01,  &
   .97035E+00, -.17586E+00, -.10184E+00,  .62649E+00, -.18503E+01,  &
   .92260E+00,  .10520E+01, -.14577E+01, -.19297E+00, -.11617E+00/

data ((ulim1(ig,ib),ig=1,3),ib=1,8)/  &
   1000.,   10000.,   1000.,  &
   10000.,  1000.,    1000.,  &
   1000.,   1000.,    2.,  &
   1000.,   1000.,    1000.,  &
   1000.,   100.,     1000.,  &
   1000.,   1000.,    1000.,  &
   1000.,   1000.,    0.3,  &
   1000.,   110.,     1000./
data (wlenlo1(ib),ib=1,8)/  &
   1.529999, .6999999, .2451000, 20.00000,  &
   12.50000, 8.333299, 9.008999, 4.642000/
data (wlenhi1(ib),ib=1,8)/  &
   4.642000, 1.529999, .6999999, 104.5149,  &
   20.00000, 9.008999, 10.30930, 8.333299/
data (solar1(ib),ib=1,3)/ 162.8590546, 556.7694702, 640.6521606/
data (ralcs1(ib),ib=1,8)/  &
  0.2306902E-09, 0.4273218E-08, 0.7142404E-07, 0.5187161E-04,  &
  0.4707304E-05, 0.8822867E-06, 0.5892014E-06, 0.4331474E-06/
data (nuum1(ib),ib=1,8)/ 0, 0, 0, 0, 0, 1, 1, 0/
data (a01(ib),ib=1,8)/  &
    .0000000E+00,   .0000000E+00,   .0000000E+00,   .2380230E+01,  &
    .3292058E+02,   .1123225E+02,   .2146121E+01,  -.8898108E+02/
data (a11(ib),ib=1,8)/  &
    .0000000E+00,   .0000000E+00,   .0000000E+00,  -.1055970E+00,  &
   -.4862469E+00,  -.7751060E-01,   .2143994E-01,   .1392107E+01/
data (a21(ib),ib=1,8)/  &
    .0000000E+00,   .0000000E+00,   .0000000E+00,   .1012857E-02,  &
    .2017528E-02,  -.3353665E-03,  -.4981068E-03,  -.7344463E-02/
data (a31(ib),ib=1,8)/  &
    .0000000E+00,   .0000000E+00,   .0000000E+00,  -.8937388E-06,  &
   -.8957686E-06,   .2563034E-05,   .1815480E-05,   .1314578E-04/

data ((npsb1(ig,ib),ig=1,3),ib=1,8)/  &
    5,    3,    0,  &
    5,    1,    0,  &
    1,    1,    5,  &
    5,    0,    0,  &
    4,    4,    0,  &
    2,    1,    0,  &
    2,    1,    3,  &
    5,    1,    0/

data ((prf1(ig,ib),ig=1,3),ib=1,8)/  &
   .1013250E+06,   .8612625E+05,   .0000000E+00,  &
   .1013250E+06,   .8612625E+05,   .0000000E+00,  &
   .1013250E+06,   .8612625E+05,   .3039750E+04,  &
   .5066250E+05,   .0000000E+00,   .0000000E+00,  &
   .8612625E+05,   .1013250E+05,   .0000000E+00,  &
   .8612625E+05,   .5066250E+05,   .0000000E+00,  &
   .8612625E+05,   .5066250E+05,   .1013250E+05,  &
   .7599375E+05,   .5066250E+05,   .0000000E+00/

data ((trf1(ig,ib),ig=1,3),ib=1,8)/  &
   .2817000E+03,   .2558000E+03,   .0000000E+00,  &
   .2817000E+03,   .2558000E+03,   .0000000E+00,  &
   .2817000E+03,   .2558000E+03,   .2299000E+03,  &
   .2558000E+03,   .0000000E+00,   .0000000E+00,  &
   .2817000E+03,   .2299000E+03,   .0000000E+00,  &
   .2817000E+03,   .2558000E+03,   .0000000E+00,  &
   .2817000E+03,   .2558000E+03,   .2040000E+03,  &
   .2817000E+03,   .2558000E+03,   .0000000E+00/

data ((xp1(ig,ik,1),ig=1,3),ik=1,7)/  &
   .1247157E-01,   .5996504E-04,   .0000000E+00,  &
   .2623545E+01,   .8041044E-02,   .0000000E+00,  &
   .8128411E-05,   .6263048E+00,   .0000000E+00,  &
   .1484021E+00,   .0000000E+00,   .0000000E+00,  &
   .7244654E-03,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((xp1(ig,ik,2),ig=1,3),ik=1,7)/  &
   .4487467E-02,   .2326860E-04,   .0000000E+00,  &
   .5147366E+00,   .0000000E+00,   .0000000E+00,  &
   .5490468E-05,   .0000000E+00,   .0000000E+00,  &
   .5409482E-03,   .0000000E+00,   .0000000E+00,  &
   .3205971E-01,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((xp1(ig,ik,3),ig=1,3),ik=1,7)/  &
   .1106090E-04,   .1005030E-04,   .1914260E+05,  &
   .0000000E+00,   .0000000E+00,   .5794290E+03,  &
   .0000000E+00,   .0000000E+00,   .7177940E+02,  &
   .0000000E+00,   .0000000E+00,   .1870680E+01,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((xp1(ig,ik,4),ig=1,3),ik=1,7)/  &
   .1790956E+00,   .0000000E+00,   .0000000E+00,  &
   .1755604E+01,   .0000000E+00,   .0000000E+00,  &
   .4080156E-02,   .0000000E+00,   .0000000E+00,  &
   .3145180E+02,   .0000000E+00,   .0000000E+00,  &
   .2784468E-01,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((xp1(ig,ik,5),ig=1,3),ik=1,7)/  &
   .5217378E-03,   .6755403E+00,   .0000000E+00,  &
   .1272084E+01,   .3978131E-01,   .0000000E+00,  &
   .5748921E-02,   .3571308E-03,   .0000000E+00,  &
   .5702444E-01,   .1850978E+02,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((xp1(ig,ik,6),ig=1,3),ik=1,7)/  &
   .1298981E-03,   .5129330E-04,   .0000000E+00,  &
   .4006586E-02,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((xp1(ig,ik,7),ig=1,3),ik=1,7)/  &
   .2897531E-04,   .5129330E-04,   .4995610E+02,  &
   .1949597E-02,   .0000000E+00,   .4341897E-03,  &
   .0000000E+00,   .0000000E+00,   .1086045E+02,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((xp1(ig,ik,8),ig=1,3),ik=1,7)/  &
   .6476305E-03,   .1053610E-03,   .0000000E+00,  &
   .6949131E+00,   .0000000E+00,   .0000000E+00,  &
   .7677455E-02,   .0000000E+00,   .0000000E+00,  &
   .8363310E+01,   .0000000E+00,   .0000000E+00,  &
   .8868708E-01,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/

data ((wght1(ig,ik,1),ig=1,3),ik=1,7)/  &
   .1415609E+00,   .7785574E+00,   .0000000E+00,  &
   .7056218E-01,   .1693635E+00,   .0000000E+00,  &
   .4375873E+00,   .5207906E-01,   .0000000E+00,  &
   .1332501E+00,   .0000000E+00,   .0000000E+00,  &
   .2170395E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((wght1(ig,ik,2),ig=1,3),ik=1,7)/  &
   .1096637E+00,   .1000000E+01,   .0000000E+00,  &
   .3161782E-01,   .0000000E+00,   .0000000E+00,  &
   .5323901E+00,   .0000000E+00,   .0000000E+00,  &
   .2418084E+00,   .0000000E+00,   .0000000E+00,  &
   .8451997E-01,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((wght1(ig,ik,3),ig=1,3),ik=1,7)/  &
   .1000000E+01,   .1000000E+01,   .5640000E-03,  &
   .0000000E+00,   .0000000E+00,   .1086900E-01,  &
   .0000000E+00,   .0000000E+00,   .1243200E-01,  &
   .0000000E+00,   .0000000E+00,   .1844170E+00,  &
   .0000000E+00,   .0000000E+00,   .7917180E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((wght1(ig,ik,4),ig=1,3),ik=1,7)/  &
   .2591315E+00,   .0000000E+00,   .0000000E+00,  &
   .3134250E+00,   .0000000E+00,   .0000000E+00,  &
   .5277705E-01,   .0000000E+00,   .0000000E+00,  &
   .1492298E+00,   .0000000E+00,   .0000000E+00,  &
   .2254367E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((wght1(ig,ik,5),ig=1,3),ik=1,7)/  &
   .3616211E+00,   .1671412E+00,   .0000000E+00,  &
   .5159410E-01,   .2359013E+00,   .0000000E+00,  &
   .3750388E+00,   .5389559E+00,   .0000000E+00,  &
   .2117461E+00,   .5800171E-01,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((wght1(ig,ik,6),ig=1,3),ik=1,7)/  &
   .8236111E+00,   .1000000E+01,   .0000000E+00,  &
   .1763889E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((wght1(ig,ik,7),ig=1,3),ik=1,7)/  &
   .8024495E+00,   .1000000E+01,   .1427190E+00,  &
   .1975506E+00,   .0000000E+00,   .5857329E+00,  &
   .0000000E+00,   .0000000E+00,   .2726962E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((wght1(ig,ik,8),ig=1,3),ik=1,7)/  &
   .1278349E+00,   .1000000E+01,   .0000000E+00,  &
   .2494413E+00,   .0000000E+00,   .0000000E+00,  &
   .2109432E+00,   .0000000E+00,   .0000000E+00,  &
   .1522931E+00,   .0000000E+00,   .0000000E+00,  &
   .2594875E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/

data ((alpha1(ig,ik,1),ig=1,3),ik=1,7)/  &
   .8242220E+00,   .9999895E+00,   .0000000E+00,  &
   .6627719E+00,   .5184072E+00,   .0000000E+00,  &
   .9999794E+00,   .6601071E+00,   .0000000E+00,  &
   .7781165E+00,   .0000000E+00,   .0000000E+00,  &
   .7519682E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((alpha1(ig,ik,2),ig=1,3),ik=1,7)/  &
   .7693452E+00,   .1385630E+00,   .0000000E+00,  &
   .5131682E+00,   .0000000E+00,   .0000000E+00,  &
   .3494183E+00,   .0000000E+00,   .0000000E+00,  &
   .7714978E+00,   .0000000E+00,   .0000000E+00,  &
   .6381509E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((alpha1(ig,ik,3),ig=1,3),ik=1,7)/  &
  -.2816690E+00,   .1000000E+01,  -.2044000E-02,  &
   .0000000E+00,   .0000000E+00,   .7768400E-01,  &
   .0000000E+00,   .0000000E+00,  -.2296670E+00,  &
   .0000000E+00,   .0000000E+00,   .9945000E-01,  &
   .0000000E+00,   .0000000E+00,   .1000000E+01,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((alpha1(ig,ik,4),ig=1,3),ik=1,7)/  &
   .9998096E+00,   .0000000E+00,   .0000000E+00,  &
   .8934920E+00,   .0000000E+00,   .0000000E+00,  &
   .4805090E+00,   .0000000E+00,   .0000000E+00,  &
   .7776389E+00,   .0000000E+00,   .0000000E+00,  &
   .7963042E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((alpha1(ig,ik,5),ig=1,3),ik=1,7)/  &
   .2766540E+00,   .7073100E+00,   .0000000E+00,  &
   .5560525E+00,   .7701662E+00,   .0000000E+00,  &
   .6394594E+00,   .2286618E+00,   .0000000E+00,  &
   .9728559E+00,   .4597031E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((alpha1(ig,ik,6),ig=1,3),ik=1,7)/  &
   .9998614E+00,   .7833650E+00,   .0000000E+00,  &
   .4105807E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((alpha1(ig,ik,7),ig=1,3),ik=1,7)/  &
   .9986984E+00,   .2954650E+00,   .2158372E+00,  &
   .3831264E+00,   .0000000E+00,   .3950504E-06,  &
   .0000000E+00,   .0000000E+00,   .4674741E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((alpha1(ig,ik,8),ig=1,3),ik=1,7)/  &
   .9913880E+00,   .5113000E-02,   .0000000E+00,  &
   .8563318E+00,   .0000000E+00,   .0000000E+00,  &
   .9979001E+00,   .0000000E+00,   .0000000E+00,  &
   .7047079E+00,   .0000000E+00,   .0000000E+00,  &
   .8926126E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/

data ((beta1(ig,ik,1),ig=1,3),ik=1,7)/  &
  -.1712392E+01,  -.3158334E+00,   .0000000E+00,  &
  -.3852223E+00,  -.5670782E+00,   .0000000E+00,  &
  -.1378367E+01,  -.7933998E+00,   .0000000E+00,  &
  -.7341006E+00,   .0000000E+00,   .0000000E+00,  &
  -.1588977E+01,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((beta1(ig,ik,2),ig=1,3),ik=1,7)/  &
  -.3969185E+00,   .4752930E+00,   .0000000E+00,  &
  -.3413553E+00,   .0000000E+00,   .0000000E+00,  &
  -.1065748E+01,   .0000000E+00,   .0000000E+00,  &
  -.5795957E+00,   .0000000E+00,   .0000000E+00,  &
  -.3413572E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((beta1(ig,ik,3),ig=1,3),ik=1,7)/  &
   .4186570E+00,  -.1006700E+01,   .4999120E+00,  &
   .0000000E+00,   .0000000E+00,   .4854630E+00,  &
   .0000000E+00,   .0000000E+00,   .4645810E+00,  &
   .0000000E+00,   .0000000E+00,  -.2546340E+00,  &
   .0000000E+00,   .0000000E+00,   .1000000E+01,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((beta1(ig,ik,4),ig=1,3),ik=1,7)/  &
  -.1718928E+01,   .0000000E+00,   .0000000E+00,  &
  -.1005695E+01,   .0000000E+00,   .0000000E+00,  &
   .9999878E+00,   .0000000E+00,   .0000000E+00,  &
  -.7752110E+00,   .0000000E+00,   .0000000E+00,  &
  -.3295007E+01,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((beta1(ig,ik,5),ig=1,3),ik=1,7)/  &
  -.2055048E+01,  -.2231248E+01,   .0000000E+00,  &
  -.4119874E+01,  -.4965271E+01,   .0000000E+00,  &
  -.2822644E+01,  -.1611550E+01,   .0000000E+00,  &
  -.3811633E+01,  -.1467305E+01,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((beta1(ig,ik,6),ig=1,3),ik=1,7)/  &
  -.6162320E+01,  -.4473330E+01,   .0000000E+00,  &
  -.2685515E+01,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((beta1(ig,ik,7),ig=1,3),ik=1,7)/  &
  -.9999836E+01,  -.5629570E+01,   .6873658E-01,  &
  -.2214210E+01,   .0000000E+00,  -.1255250E+01,  &
   .0000000E+00,   .0000000E+00,  -.1283144E+01,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/
data ((beta1(ig,ik,8),ig=1,3),ik=1,7)/  &
  -.4327339E-01,   .4952220E+00,   .0000000E+00,  &
  -.2665102E+01,   .0000000E+00,   .0000000E+00,  &
  -.6325287E+01,   .0000000E+00,   .0000000E+00,  &
  -.1731381E+01,   .0000000E+00,   .0000000E+00,  &
  -.4994114E+01,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00,  &
   .0000000E+00,   .0000000E+00,   .0000000E+00/

solc = 1370.
nb = 8
ng = 3
ns = 3
ngb = 2

do ip = 1,13
   icat = kkat(ip)
   ignu = max(1,min(2,nint(gnu(icat))))
   do ib = 1,8
      do ic1 = 1,5
         ocoef(ic1,ib,ip) = ocoef1(ic1,ib,ip,ignu)
      enddo
      do ic2 = 1,2
         bcoef(ic2,ib,ip) = bcoef1(ic2,ib,ip,ignu)
      enddo
   enddo
enddo

do ip = 1,7
   ignu = max(1,min(2,nint(gnu(ip))))
   do ib = 1,8
      do ic3 = 1,5
         gcoef(ic3,ib,ip) = gcoef1(ic3,ib,ip,ignu)
      enddo
   enddo
enddo

do ib = 1,8
   do ig = 1,3
     ulim(ig,ib) = ulim1(ig,ib)
     prf(ig,ib) = prf1(ig,ib)
     trf(ig,ib) = trf1(ig,ib)
     npsb(ig,ib) = npsb1(ig,ib)
   enddo
enddo

do i = 1,nb
   wlenlo(i) = wlenlo1(i)
   wlenhi(i) = wlenhi1(i)
   ralcs(i) = ralcs1(i)
   nuum(i) = nuum1(i)
   a0(i) = a01(i)
   a1(i) = a11(i)
   a2(i) = a21(i)
   a3(i) = a31(i)
enddo

do is = 1,ns
   solar(is) = solar1(is)
enddo

do ib = 1,8
   do ig = 1,3
      do ik = 1,7
         xp(ig,ik,ib) = xp1(ig,ik,ib)
         wght(ig,ik,ib) = wght1(ig,ik,ib)
         alpha(ig,ik,ib) = alpha1(ig,ik,ib)
         beta(ig,ik,ib) = beta1(ig,ik,ib)
      enddo
   enddo
enddo

!     do the solar, rayleigh scatter, and continuum abs.
do ib=1,nb

   if (ib.le.ns) then
      solar(ib) = sunavg(1.e4/wlenhi(ib),1.e4/wlenlo(ib),solc)
      call rayleigh(wlenlo(ib),wlenhi(ib),ralcs(ib))
   else
      call csband(wlenlo(ib),wlenhi(ib),ralcs(ib))
   endif

enddo

return
end subroutine harr_radinit1
